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ABSTRACT 

The aim of the present work is to show that, contrary to popular belief, galaxy clusters 
are not expected to be self-similar, even when the only energy sources available are 
gravity and shock-wave heating. In particular, we investigate the scaling relations 
between mass, luminosity and temperature of galaxy groups and clusters in the absence 
of radiative processes. Theoretical expectations are derived from a polytropic model 
of the intracluster medium and compared with the results of high-resolution adiabatic 
gasdynamical simulations. It is shown that, in addition to the well-known relation 
between the mass and concentration of the dark matter halo, the effective polytropic 
index of the gas also varies systematically with cluster mass, and therefore neither the 
dark matter nor the gas profiles are exactly self-similar. It is remarkable, though, that 
the effects of concentration and polytropic index tend to cancel each other, leading to 
scaling relations whose logarithmic slopes roughly match the predictions of the most 
basic self-similar models. We provide a phenomenological fit to the relation between 
polytropic index and concentration, as well as a self-consistent scheme to derive the 
non-linear scaling relations expected for any cosmology and the best-fit normalizations 
of the M-T, L-T and F-T relations appropriate for a ACDM universe. The predicted 
scaling relations reproduce observational data reasonably well for massive clusters, 
where the effects of cooling and star formation are expected to play a minor role. 
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1 INTRODUCTION 

The physics of massive galaxy clusters is relatively simple, 
at least compared to that of smaller objects. In the stan- 
dard cold dark matter (CDM) scenario, their mass is domi- 
nated by the dark component, while most baryons are in the 
form of a hot diffuse plasma in hydrostatic equilibrium with 
the gravitational potential created by the CDM halo. The 
intracluster medium (ICM) gas is shock-heated to approxi- 
mately the virial temperature of the object, and its thermal 
bremsstrahlung emission has been detected by X-ray satel- 
lites for the last three decades. In the absence of any other 
process, it is often stated that galaxy clusters are expected 
to be self-similar, and their global properties should obey 
power-law scaling relations. 

As long as the shape of a cluster's potential does not 
depend systematically on its mass, the radial structure of 
the ICM ought to be scale-free, and the global properties 
of galaxy clusters, such as halo mass, emission- weighted 
temperature, or X-ray luminosity, would scale self-similarly 
lKaiserl ll98fih. Indeed, numerical simulations that include 
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adiabatic gasdynamics are reported to prod uce clusters of 



galaxies that obey such scaling laws (e.g. | Na varro et al 
|l99||lEvrard et alJll996t iBrvan fc NormarJll998tlEke et a* 

In real clusters, deviations fro m self-similarity ar e ex- 
pected to arise from merging (e.g. Ijing fc Suto 2000J) and 
additional physics act i ng on the intracluster gas (see e.g . 
iTozzi fc NormarJl200ll: iBabul et all 12001 IVoit et al] 120021 
and references therein). Radiative cooling and energy injec- 
tion by supernova and/or active galactic nuclei (AGN) may 
be particularly relevant for low-mass systems, where they 
can make a significant contribution to the total energy bud- 
get. 

Observations seem to corroborate that the self-similar 
picture is indeed too simplistic, and it fails to predict the 
observed scalings of cluster mass and luminosity with respect 
to the ICM gas temperature. 

Although some observational studies (e.g. [Ho rner et al] 
ll999HNeumann fc Arnaudlll999l:lAllen et al]l200ll) are c~ 
sistent with the self-similar expectation, M oc T 1 ' 5 , the ob- 
served mass-te mperature relation has often been found to 
be steeper (e.g. ISanderso n et al.l l200ot). particularly in the 
group regime (e.g lNevalamen^taT . 2000; Finoeuenov et al] 
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120011: IXu et~ai]l200ll; lArnaud et ai]l2005l) . It has also been 
noted fcttori et al.ll2002l) that the slope of the M — T re- 
lation might depend on the limiting overdensity. Regarding 
the normalization, the reported value is in most cases ~ 40 
per cent lower than in numerical simulations, although there 
is a certain degeneracy with the precise value of the slope. 

On the other hand, it has been k nown since the 
first generation of X-r ay satellites (e.g. lEdee fc Stewart! 
119911: iDavid et al.lll993ft that the slope of the luminosity- 
temperature relation, a ~ 3, is significantly steeper 
than the self-similar value, a = 2 (although the ex- 
act value v aries for a limited en e rgy band). It has been 
shown (e.g. lAllen fc Fa,bianlll99st iMarkevitch et all Il998l 
lArnaud fc Evrardl ll999) that the scatter in the L — T re- 
lation is significantly reduced, and the discrepancy is some- 
what less severe, when cooling flows are excised or samples 
with only weak cooling cores are considered . Actually, it has 
been recently claimed dO'Hara et al.l l200fJ) that cool core 
related phenomena, and not merging processes, are the pri- 
mary contributor to the scatter in all the scaling relations. 

Observational data have thus motivated significant ef- 
forts attempting to build a physical model of the ICM 
that breaks self-similarity, either by removing low-entropy 
gas from th e centres of clusters via radiative cooling 
1 Bryan| ^ QQOj) or bv introducing no n-gravitational heating 
l Evrard fc Henrvll99]]:lKaiserlll99J) . In both cases, the 'ex- 
cess' entropy produces a flattening of the density profile that 
brings the X-ray properties of the modelled cluste rs in agree- 
ment with the observed s caling relations (see e.g. IVoit et ail 
kOO&lBorganT et al.l2004l and references therein) . Neverthe- 
less, the source and precise amount of heating and cooling 
requ ired are still a matt er of debate. Recent observations 
(e.g. iPonman et al]l2003l) suggest that the shape of the en- 
tropy profile is similar in groups and clusters of galaxies, 
which rules out the simplest scenarios. 

In this paper, we claim that dark matter haloes are not 
exactly self-similar, and therefore both the cluster's poten- 
tial and the properties of the ICM gas do indeed depend on 
the total mass (or temperature) of the object. In particular, 
there is no compelling reason to expect that the scaling re- 
lations between any two physical properties, integrated up 
to a given overdensity, should obey a power law, even in the 
purely gravitational case. 

We present a theoretical prediction of these relations 
based on a polytropic model of the ICM, and compare it 
with a set of high-resolution adiabatic gasdynamical simu- 
lations. It will be shown that self-similar models implicitly 
assume that all clusters have the same concentration and 
polytropic index. Relaxing these hypotheses yields the scal- 
ing relations derived in Section |5] which we compare with 
the results of numerical experiments in Sectional Observa- 
tional implications are discussed in Sectional and Section |S] 
summarizes our main conclusions. 



2 THEORETICAL MODEL 



As shown by lAscasibar et alJ <l2003l) , relaxed clusters and 
minor mergers found in adiabatic gasdynamical simulations 
can be considered to be in approximate thermally-supported 
hydrostatic equilibrium up to ~ 0.8r2oo- Furthermore, the 
ICM gas is fairly well described by a polytropic equation 



of state with an effective polytropic index 7 ~ 1.18. Usingj 
the p henomenological formula proposed by iNavarro et alJ 
( 1997, hereafter NFW) to model the density profile of the 
dark matter halo, 



p(r) = 



(r/r s )(l + r/r s )2' 
the gas temperature is given by 
ln(l + r/r s ) 



T(r) = To 



where the central temperature, 



fcTb = 47rG/im p 



7" 



7 



1 2 
-p s r s , 



(1) 



(2) 



(•3) 



is set by the boundary condition that the ICM density and 
temperature vanish at infinity. The gas density profile can 
be computed from the polytropic relation 



Ps(r) = Po(7) 



ln(l + r/rs) 
r/r s 



(4) 



where the central gas density, po(7), can be constrained by 
normalizing the baryon fraction to match the cosmic value 
at x h = r h /r s ~ 3, 



f \ Qb 
Po[V = 75 — Ps 



ln(l +T b ) 
x b 



(5) 



T his analytic prescriptio n is simpler than the one proposed 
in lAscasibar et all £2003), and it provides better results for 
extreme values of the polytropic index 7. The choice x b ~ 3 
is somewhat arbitrary, and in principle one could leave the 
normalization of the gas density as a free parameter of the 
model. 

We consider, though, that it is desirable to reduce the 
number of free parameters as much as possible. In fact, the 
very existence of relatively tight scaling relations suggests 
that real galaxy clusters can indeed be described by only 
one free parameter, which could be taken to be the mass 
of the halo. Once Xb is set, our model still has three pa- 
rameters, two of them related to the dark matter halo (the 
characteristic density and radius, p s and r s ) and one re- 
lated to the intracluster gas (the effective p olytropic index 
7). The first two are known to be correlat ed jNavarro et alJ 
Il997t iBullock et aljEobH: lEke et alJl200ll) . and we propose 
a phenomenological relation between polytropic index and 
concentration in Section [3] below. 

It is important to note, though, that the fact that clus- 
ters could be described by a one-parameter family of func- 
tions would give rise to universal scaling relations both for 
their radial profiles as well as for their global (integrated 
or averaged) physical properties. However, it does not im- 
ply that such relations ought to be self-similar in any sense. 
The precise functional form of the different scalings would 
be specified by the two independent relations between p s , r s 
and 7. 

In the present work, we are interested in the scaling 
relations between several quantities, integrated up to the 
radius ta encompassing an overdensity A with respect to 
the critical density, i.e. 



47T -i 

Ma = M(r A ) = A "g-Pc r A , 



(6) 
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2.8 x 10 11 h 2 M© Mpc -3 is the critical 



where p c = ^§ 
density and Ho = 100 h km s Mpc 1 is the Hubble con- 
stant. Defining ca = ta/t s and g(x) = [ln(l + x) — jtj] , 
we obtain 



47rp s r 3 g 1 (c A ) 



M A 
and 

AHn 3 , » 

Ps = T"7=r c A s(ca). 

07TLr 

Let us define the function 



f 


ln(l + a;)" 


Jo 


X 



x 2 dx, 



(7) 



(8) 



(9) 



which must be integrated numerically. In terms of this func- 
tion, we can express in a very compact form both the cumu- 
lative gas mass, 



Analogously, expressions Q and 1101 tell us that the 
cumulative baryon fraction does not depend explicitly on 
the object mass or temperature, 



Fa 



Ma n b 



po (7) Q d E 
fib 



5(ca) y(- 



1 



, CA) 



(18) 



p s I ib 7 
Finally, the luminosity-temperature relation can be ex- 
pressed in the form 

r A Ax-ffo 



( / um p )2 2'i ttG 2 



(kTi) y L x(7,c A ) 



with 



iLx(7, ca) 



fxtrr) 



7 ca 



S 2 (ca) 



y 3 (^r + jcA) 
y 2 (^T + |,CA) 



(19) 



(20) 



according to JSJ, ltT4l and ^15t . 



M g A = 47rp r s 3 t/( 



,CA) 



7-1' 

and mass-weighted temperature, 



Ta=T 



y(rrr,c A ) 



y( 



7-1 ' 



ca) 



(10) 



(11) 



Assuming that thermal bremsstrahlung is the dominant 
cooling mechanism, the X-ray power radiated by the ICM 
gas per unit volume may be estimated as 



647T /7T\ 



V— T 

\4vreo J 



kT 



(m e c 2 ) 3 



gn, 



J2 z * n *> ( l2 ) 



where e, m B and n e are the electron charge, mass and num- 
ber density, respectively, eo is the permittivity of free space, 
k is Boltzmann's constant, c is the speed of light and g is 
the average Gaunt factor, which we take to be unity. The 
sum takes into account the atomic number Zi and number 
density m of each ion species i. For a fully-ionized plasma 
of primordial composition (~ 75 per cent of the mass in 
hydrogen and 25 per cent in helium), 

( 1 Vf— Y 

lM Q Mpc" 3 / VkeVy 



2 x 10 



erg s 



"'Mpc -3 . (13) 



Integrating up to r a , the bolometric X-ray luminosity would 
be 



£x = Ax pi (fcT )5 47rr s 3 y(- 



2 1 
— +2,c A ) 



with A x ~ 2 x 10 17 erg s" 1 Mpc 3 M " 2 keV" 
emission-weighted temperature is given by 

2/(^T + |,ca) 



t£ = t 



2/(^t + g, ca)' 



(14) 

f , while the 
(15) 



Combining equations 0, JHJ and I15L simple al- 
gebra yields the mass-temperature relation 



M A 



V2 
GH 



-1/2 



prrip 



3/2 



yMT(7,CA), 



where 



Vmt(7j ca) 



7-1 , M^=I + 

CA g(CA) 



CA) 



y(^h+ 2,CA) 



(16) 



(17) 



3 SIMULATIONS 

For constant values of the polytropic index and concentra- 
tion, equations 1161 . 11811 and 119H become the well-known 
self-similar scalings, with logarithmic slopes 3/2, and 2, re- 
spectively. The precise 'universal' values of 7 and ca would 
simply set the normalization. 

However, both quantities might well depend systemati- 
cally on the mass of the object. We address such dependence 
in the present section, where we also compute the expected 
scaling relations and compare them to our numerical data. 



3.1 Numerical experiments 

Our cluster sample consists of 42 objects formed in a 
flat ACDM universe (fi m = 0.3; Q h = 0.04; Q A = 0.7; 
h = 0.7; us = 0.9). 28 of them have been extracted 
from a 80 h" 1 Mpc cubic box simula ted with a version o f 
the parallel Tree-SPH code Gadget JSpringel et al.ll200ll) 
that implements the entropy -conserving scheme proposed by 
ISpringel fc Hernauistl i2002T) . For a thoro ugh description o f 
these experiments, the reader is referred to lAscasibari (2003). 
In order to extend our numerical sample of clusters to a 
wider temperature (mass) range, we have also simulated a 
500 h' 1 Mpc box with the code Gadget2 ISpringelll2o"o"5T) . 
from which we have considered 14 objects. Details about 
these simulations can be found in lYepes et alj j2004l) . 

In each case, high resolution has been achieved by 
mean s of the multiple-mass technique (see iKlypin et alJ 
2001). An unconstrained random realization of the ACDM 
power spectrum was generated with 1024 3 and 2048 3 par- 
ticles for the 80 and 500 Mpc boxes, respectively. 
Haloes were selected at z — from a low-resolution ex- 
periment evolved with 128 3 dark matter particles, and then 
re-simulated with three and five levels of mass refinement 
(so that the final mass resolution of both subsamples is 
similar). The gravitational softening length was set to e — 
2 — 5 h^ 1 kpc, depending on number of dark ma tter particles 
withi n the virial radius of the object f following iPower et alJ 
2003). Gas particles have only been added in the highest 
refinement level. 

Basic information about the objects in our numerical 
cluster sample is summarized in Tabled They span two or- 
ders of magnitude in mass (~ 10 13 — 10 15 Mq) and cover 
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Table 1. Description of our cluster sample. Number of gas particles within r200, physical properties at overdensities A = 2500, 500 and 
200, best-fitting characteristic density (in units of the critical density), radius (in h _1 kpc) and effective polytropic index. Masses are 
expressed in 10 13 Mq, temperatures in keV and X-ray luminosities in 10 44 h erg s _1 . The baryon fraction is given in units of the cosmic 
value, Q b /f! m . 
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10.11 

2.14 
1.62 
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10.64 
8.14 
6.14 
5.74 
10.99 
3.29 
4.40 
7.38 
1.43 
1.92 
1.17 
1.32 
0.99 
0.68 
0.69 
0.98 
1.51 
0.64 
0.43 
1.88 
0.81 
3.46 
1.06 
0.83 
1.20 
1.21 
0.86 
0.79 
4.23 
4.24 
1.28 
2.66 
2.75 
3.07 
2.83 
0.80 
1.72 
1.55 



a temperature range between 0.5 and 11 keV. The total 
number of gas particles within the virial radius is always 
iVg as ^ 2 x 10 4 , the number of dark matter particles be- 
ing slightly higher (inversely proportional to the cumulative 
baryon fraction, F). 

For each object, the centre of mass was found by an 
iterative procedure. Starting with an initial guess, we com- 
pute the centre of mass within a sphere of 500 hT 1 kpc. 
The sphere is moved to the new centre until convergence is 
reached. The radius of the sphere is then decreased by 10 
per cent, and the process continues until the sphere contains 
200 dark matter particles. The radii r2soo, fsoo and r2oo are 
obtained from the overdensity profile around the final centre 
of mass. The total mass, baryon fraction, X-ray luminosity 



and emission- weighted temperature quoted in Tabled have 
been computed within those radii. 



3.2 Concentration and polytropic index 

The values of the parameters ca and 7 have been com- 
puted by means of a global fit to the gas density, dark mat- 
ter density, total mass and effective polytropic index (i.e. 
gas temperature versus gas density) profiles, averaged over 
logarithmically-spaced spherical shells between 0.1r2oo and 
?"200- More precisely, we minimize the quantity 

X = X P + Xm + Xt + X 2 ~t , (21) 
where 
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m 20 o [ h ~ 1 M ] 

Figure 1. Concentration at different overdensities and ef- 
fective polytropic index of each object. Solid lines represent 
our th eoretical model, bas ed on the mass-concentration rela- 
tion of iBullock et all 1200 lfl and our proposed fit to the depen- 
dence of the effective polytropic index on concentration, equa- 
tion 1231 , Dotted lines show the one -sigma scatter expected from 
Ac vir /c vir ~ 0.3 ICoh'n et alj|2004l). Th e linea r fit to the 7 - c 
relation proposed bv TKoniateu^^eTialJ 1200 ll) is shown by the 
dashed lines on the right panels. 



"bins 

E 

6=1 



log 2 [x(b)/a: modc i(b)] 
l/n(b) 



(22) 



and x denotes the number of particles within each bin, 
n(b), the total enclosed mass, M(b), the average tempera- 
ture within the bin, T(b) = X^T-i ^V^Wi an( ^ the quantity 
T(6)n 1 " 7 (6). 

A grid of analytical profiles is generated for the intervals 
1.1 < 7 < 1.25 and 10 < r s /(/i _1 kpc) < 1000, in uniform 
steps A7 = 0.001 and Ar s — 1 h~ kpc. Since all the mea- 
sured quantities are proportional to the characteristic den- 
sity p a , its best-fitting value has been trivially found from 
the average of the logarithmic residuals. 

Results of the minimization procedure are quoted in 
Table Q and best-fitting concentrations and effective poly- 
tropic indices are plotted in Figure Solid lines depict the 
toy model for the mass-concentration relation proposed by 



IBullock et aflfcOOlh 1 . with F = 0.001 and K = 3. We trans- 
form the values of c v j r ~ C100 to the other overdensities ac- 
cording to the NFW profile. Dotted li nes show the one- sigma 
scatter Ac v ir/c v ir — 0.3 reported bv lColm et al.l ll2004) for 
relaxed systems. Finally, we find that the phenomenological 
relation 



7 = a + b C200 



(23) 



with a = 1.145 ±0.007 and b = 0.005 ±0.002 fits reasonably 
well our results for the polytropic index (although there ex- 
ists a certain degeneracy between the best-fitting values of 
both parameters). Dotted lines show the one-sigma scatter 
in 7(c) expected from Ac v i r - 

A correlation between 7 and c is expected if both 
quantities (and therefore the radial structure of the ICM) 
vary smoothly with the mass of the object. In fact, an ap- 
pr oximately linear dependen ce has already been advocated 
by iKomatsu fc Seliakl \2QQ ll) in order to enforce constant 
baryon fraction at large radii. Their fit (dashed lines in Fig- 
ure is however significantly steeper than equation I23H . 
Although it works somewhat better for large c, it does not 
seem to adequately describe the least concentrated systems, 
suggesting that, most probably, the precise functional form 
of 7(c) is not as simple as a straight line. We therefore advise 
against extrapolating our fit towards values of the concen- 
tration parameter outside the range covered by the present 
work. Moreover, additional physics is expected to play an 
important role in less massive (more concentrated) systems, 
and therefore the polytropic approximation itself will no 
longer to be valid, since it fails to describe the presence of a 
central cool core. 

Some of our objects deviate appreciably from both the 
mass-concentration relation and the 7 — c relation given by 
expression 12311 . These tend to be merging systems, which 
have formed (or are forming) more recently than relaxed ob- 
jects. In the spherical collapse picture, that means they have 
collapsed around density peaks on larger scales, and there- 
fore the resulting density profiles are less concentrated than 
relaxed h aloes of the sam e mass llAscasibar et al 1 I2004L As 
noted bv lGottlober et al] \2QQ ll) . merging is more common 
on the scale of galaxy groups. The effective polytropic index 
seems to be systematically higher in these objects, although 
it is important to bear in mind that they are not particularly 
well described by a polytropic equation of state. Actually, 
the gas distribution shows obvious asymmetries, as well as 
an offset between the gas and dark matter peak which re- 
sults in an artificially flat gas density profile in the central 
regions. Such flattening is responsible for both the abnor- 
mally low baryon fraction measured at A = 2500 and the 
unusually high value of 7 obtained by our fitting routine. 



3.3 Scaling relations 

The scaling relations of total mass, baryon fraction and X- 
ray luminosity with respect to the emission-weighted tem- 
perature are represented in Figure |2] divided by the ap- 
propriate values of the structure factors Ymt, Fa and Ylx 



1 Nearl y identical result s are obtained when the prescription 
given in lEke et alj 1200 J) is used. 
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Figure 2. M — Tx, F — Tx and Lx — Tx scaling relations of our cluster sample at different overdensities, corrected by the factors Ymt, 
Fa and Ylx, respectively. Dashed lines represent the theoretical predictions given by expressions 1161 , 1181 and 1191 , 



corresponding to each object. These are computed by sub- 
stituting the best-fitting values of 7 and ca into expressions 
O, GHJ and CT . 



As shown bv I Ascasibar et ail <l200ot) . poly tropic models 
provide a fairly accurate description of the radial structure of 
galaxy groups and clusters. It is therefore not surprising that 
they are able to match the scaling relations as well. When 
the factors Ymt, Fa and Ylx are taken into account, both 
the normalization and the logarithmic slope (equal to the 
self-similar models) of the scaling relations are correctly pre- 
dicted by equations 1161 . 1181 and 11911 . The scatter around 
the theoretical expectation is quite low, and only merging 
systems deviate appreciably from the predicted relation. In 
these objects, the dark matter potential may differ consid- 
erably from the NFW form, and the assumptions of hydro- 



static equilibrium and a poly tropic equation of stat e provide 
rather poor approximations iAscasibar et al.|l200ot) . 

Uncorrected scaling relations are plotted in Figure [3] 
Solid lines show our theore tical prediction ; usin g the mass- 
concentration relation from lBullock et al] <l200ll) and our fit 
1231 1 to estimate ca and 7 as a function of the cluster mass. It 
turns out that the dependencies on concentration and poly- 
tropic index seem to cancel each other so that the resulting 
scaling relations look as if the clusters were actually self- 
similar. Indeed, all the relations considered in the present 
study can be accurately fit by an 'average' concentration 
Cvir = 8 (which implies C200 — 6, c 5 oo — 4, C2500 — 1-8 and 
7 ~ 1.176). The corresponding scaling relations, 
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Figure 3. Numerical scaling relations (points), compared to our theoretical prediction (solid lines) based on a variable c(M) and 7(c), 
see text. Dotted lines show the scatter arising from Ac v i r /c v ; r ~ 0.3, and dashed lines indicate the scaling relations expected for a fixed 
concentration c v i r = 8 and effective polytropic index 7 = f .176. 



Table 2. Normalization of the approximate scaling relations 
obtained for c v ; r = 8 and 7 = l.f76. is expressed in 

f0 14 h- 1 Mq, L a in 10 43 h erg s~ 4 and F A in units of the 
cosmic baryon fraction. 



A 


M A 


F A 


T A 


2500 


5.65 


0.771 


2.18 


500 


5.47 


0.885 


2.69 


200 


4.71 


0.885 


2.74 



have been plotted as dashed lines, and their normalizations 
are given in Table [5] 



It is to some extent remarkable that concentration and 
polytropic index conspire to produce scaling relations that 
match so closely the self-similar slope. Such effect is illus- 
trated more clearly in Figure 0] where the contributions of 
7 and c to the scaling relations are plotted separately. 

The value of the concentration sets the ratio between 
the radii ta and the characteristic radius r s . Since massive 
objects are substantially less concentrated than smaller sys- 
tems, their r& are much closer to the centre in terms of r s , 
and thus the mass within a given overdensity (which is an 
increasing function of r/r s ) will be smaller than indicated by 
the 'average' scaling relation based on a higher value of c v ir- 
Conversely, Ma in the least massive objects would be biased 
high with respect to the self-similar relation. The effect is 
particularly noticeable for A = 2500, where ta <C r s and 
the enclosed mass is a rapidly increasing function of r/r s ; 
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Figure 4. Effects of concentration (dotted) and polytropic index (dashed) on initially self-similar scaling relations where c v 
7 = 1.176 for all objects (grey solid). Black solid lines show our prediction, taking into account both effects simultaneously. 



8 and 



for r > r s , it increases only logarithmically, and differences 
in ta/Vs are much less important. 

The emission-weighted temperature decreases with 
r/r s , and therefore it is expected to be biased high (low) for 
large (small) systems. However, the effect is arguably small, 
since the average is biased towards the central part, where 
the gas temperature is roughly constant. As shown by the 
dotted lines on the top panels of Figure |1J the net result is a 
shallower mass-temperature relation at A = 2500, while no 
significant change can be appreciated at lower overdensities. 

The baryon fraction and the X-ray luminosity are steep 
functions of r/r$, and hence more sensitive to the precise 
value of Cvir- The effect on these quantities is also stronger 
at large overdensities, but unlike the M — T relation, the sys- 
tematic variation of concentration would yield a noticeable 
imprint on the F — T and L — T relations at A = 200. 



On the other hand, the effective polytropic index con- 
trols the slope of the gas density profile. According to 123H . 
larger masses, which imply lower concentrations, also mean 
lower effective polytropic index. The gas density profile be- 
comes increasingly steep, which boosts the central baryon 
fraction and the X-ray luminosity, increasing the emission- 
weighted temperature only slightly and leaving the the to- 
tal mass mostly unaffected. Given the small variation of the 
polytropic index throughout the interesting mass range, the 
effect is never larger than a factor of two, comparable to or 
smaller than that of the concentration, but it always acts in 
the opposite sense. 

Finally, we would also like to note that, apart from the 
normalization, our model can estimate the scatter around 
the average scaling relations by assuming that it arises com- 
pletely from the scatter in the mass-concentration relation. 
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Table 3. Normalization (in 10 14 h _1 M ) and logarithmic 
slope a of the M — T relation reported in previous numerical 
studies. 



Reference 


A 


< 


a 


Navarro et al. ( 19951 


200 


6.14 


1.5 


Evrard et al. ( 1996) 


2500 


7.85 


1.5 




500 


7.84 


1.5 


Brvan & Norman (19981 


200 


8.68 


1.5 


Pen (1998) 


200 


7.28 


1.5 


Eke et al. (19981 


100 


6.12 


1.5 


Yoshikawa et al. ("20001 


100 


5.56 


1.5 


Mathiesen & Evrard ("20011 


500 


8.18 


1.52 




200 


7.65 


1.54 


Muanwone et al. (2002) 


200 


13.9 


1.5 



Such estimate, indicated by the dotted lines in Figure^ has 
been obtained by combining the value c' = 1.3c v i r with the 
effective polytropic index 7' = 7(c v ir/1.3), and c" = c v i r /1.3 
with 7" = 7(1.3c v ir), where in both cases 7(c) has been com- 
puted according to equation 12311 . Comparing to Figure[5] it 
seems that the internal structure of clusters (accounted for 
by the factors Ymt, Fa and YLx) is indeed responsible for a 
significant fraction of the scatte r in the observed sca ling re- 
lations, as recently suggested bv lO'Hara et alJ ll2005l) . Note, 
however, that in our case the differences in internal struc- 
ture are obviously not related to the presence of a cool core 
or the action of any external source of energy, but rather to 
the different formation histories of each object. 

3.4 Comparison with previous work 

The mass-temperature relation in the absence of radiative 
processes has been extensively studied by means of cosmo- 
logical numerical simulations. A summary of previous results 
is given in Tabled] Slopes consistent with 3/2 are found in 
most studies, with normalizations showing a relatively low 
scatter around Mfr ~ (7-8) x 10 14 h' 1 M©. Similar results 
are obtain ed when cooling and stellar feedback are consid- 
ered (e.g. iBorgani et al] 1200411 . although there is a trend 
towards steeper slopes and lower normalizations, in better 
agreement with observational data. 

The M — T relation predicted by our model (see Ta- 
ble |2l is also considerably lower than the values found in 
previous experiments based on purely adiabatic gasdynam- 
ics. We think that this is to a great extent a resolution effect, 
coupled to the use of an entropy-conserving scheme to solve 
the SPH equations. 

Most if not all of the earlier work on the adiabatic scal- 
ing relations of ga laxy clusters relied on the traditional for - 
mulation of SPH dLucvlll977l : iGingold fc MonagharJll977ft . 
It has been recently shown (e.g. |Springel fc Hernauist 2002; 
lAscasibar et al.ll2003l : lO'Shea et all200SH that poor entropy 
conservation leads to spurious entropy losses in the cluster 
cores, and thus previous codes tend to systematically over- 
estimate the central density and underestimate the central 
gas temperature. 

On the other hand, lack of resolution results in artifi- 
cially flattened density and temperature profiles, i.e. nearly 
isothermal and isentropic cores. Actually, some stu dies (e.g. 
iMathiesen fc EvrardlEoOlT: iMuanwong et al.ll2002l) even re- 



Table 4. L — T relation found in previous simulations, with 
expressed in 10 43 h erg s — 1 . 



Reference 


A 


T A 

L o 


a 


Navarro et al. ("19951 


200 


4.61 


2 


Brvan & Norman (19981 


200 


2.67 


2 


Eke et al. ( 19981 


100 


0.85 


2 


Bialek et al. (20011 


500 


2.44 


2.02 



port decreasing temperature profiles towards the centre. 
This is in strong disagreement with our results, as well as 
with those of independent numerica l work based on h igh- 
resolution Eulerian simulations (e.g. lLoken et al1l2002t) . in 
which the temperature profile in the absence of radiative 
processes is also found to decrease monotonically with ra- 
dius. Under such conditions, the emission-weighted temper- 
ature (which is biased towards the central, dense and X- 
ray bright regions of the cluster) is larger than the mass- 
weighted average, and thus the resulting normalization of 
M — T relation becomes considerably lower when Tx is used 
and the central parts of the objects under study are well 
resolved. 

Concerning the Lx — Tx relation, there is relatively lit- 
tle numerical work based on adiabatic gasdynamical simu- 
lations. This is also a reflection of the stringent resolution 
requirements, due to the fact that a significant fraction of 
the X-ray photons are expected to be produced in the in- 
nermost regions. A further problem affecting cosmological 
numerical experiments is that the smallest objects are typ- 
ically resolved with less particles, so their bolometric X-ray 
luminosity is underestimated and t he resulting L — T rela- 
tion is artificially steepened (see e.g. [Brro, n fc Normanll998l : 
lYoshikawa et al J 120001: lYepes et alJl2004D . Most of our ob- 
jects (see Table have more than 10 5 gas particles, and in 
principle the y should not be severely af fected by this prob- 
lem (see e.g. IBorgani et aljEooj . I2005H . Nevertheless, it is 
always wise to bear this consideration in mind when drawing 
conclusions from numerical data. 

Table|3]shows several fits to the adiabatic L — T relation 
reported in the literature. The normalizations are in this 
case broadly consistent with our results (Table although 
the scatter between different estimates is extremely large. 
This is not entirely unexpected, given the sensitivity of the 
X-ray luminosity to the details of the gas density profile in 
the central regions. 

Finally, the baryonic content of galax y cluste rs has been 
recently investigated by iKravtsov et alJ d200ot) . For their 
adiabatic simulations, they find -F2500 = 0.85 ± 0.08 and 
-F500 = 0.94 ± 0.03 in units of the cosmic value. As noted by 
these authors, the baryon fractions obt ained with the Eu- 
lerian code ART (IKravtsov et alJ 2002) arc systematically 
less concentrated than those obtained with Gadget, even 
when entropy conservation is enforced, but nevertheless the 
cumulative baryon fractions beyond T2500 a re about about 
3 — 5 per cent higher. Since the results of IKravtsov et al] 
(2005) are based on a subset of the cluster sample studied 
here, the interested reader is referred to that paper for an 
extensive comparison between both codes. 

Our res ults are also compat ible with the baryon fraction 
measured bv lEttori et alJ ll2006J) in their non-radiative runs, 
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when the standard implementation of the artificial SPH vis- 
cosity is used. As pointed out by these authors, it is interest- 
ing how an improved scheme may lower the baryon fraction 
in the innermost regions by about 15 per cent. Although 
not dramatic, it seems clear that the details of the numeri- 
cal technique do have a measurable impact on the predicted 
radial profiles near the centre, as well as on the scaling re- 
lations at high overdensities. These relatively small discrep- 
ancies between different algorithms should nevertheless be 
considered as part of the theoretical uncertainty. 



4 OBSERVATIONS 

It is the aim of the present work to provide a sound theo- 
retical prediction of the cluster scaling relations when only 
gravity, adiabatic gasdynamics and shock-wave heating act 
on the intracluster medium. In real life, additional physi- 
cal processes, such as radiative cooling, energy injection by 
stars and AGN or thermal conduction, may play an impor- 
tant role or even determine the exact form of the scaling rela- 
tions. However, the influence of all these phenomena outside 
the central regions is expected to be relatively small when 
compared to shock heating, especially for the most massive 
systems. 

Therefore, one may expect a priori that the adiabatic 
scaling relations roughly match the observed ones. Depar- 
tures would measure the effect of additional physics, and 
in principle should be more noticeable at high overdensities 
and for low-temperature systems. 

Our predictions for the adiabatic case are compared 
with observational data in Figure |S] We find a fairly good 
agreement, both in shape and, to some extent, scatter, with 
the observed M — T relation at differen t overdensities. Only 
the dataset from Piffarct ti et alJ (|2005j) is not well described 
by our model. From visual inspection of Figure though, 
it seems that the dark matter masses inferred for some of 
these systems are lower not only than our results but also 
than the other observations. 

On the other hand, the luminosity-temperature relation 
observed for galaxy clusters is roughly consistent with our 
theoretical prediction, but the slope is appreciably steeper. 
As one approaches the group regime, real systems can be 
one order of magnitude less bright than the model. The X- 
ray luminosity is much more sensitive to the details of the 
central parts than the total mass or the emission-weighted 
temperature, and thus it is not surprising that the L — T 
relation deviates from the adiabatic prediction more signifi- 
cantly than the mass-temperature relation. 

The lower X-ray emissivity seems to be intimately con- 
nected to the shape of the gas density profile. Our models 
correctly predict that the baryon fraction should be an in- 
creasing function of radius, and such a trend is clearly consis- 
tent with the observational data. However, they also predict 
that, at a given overdensity, the baryon fraction should be 
roughly independent on cluster mass. 

This is blatan t ly at odds with observations. As noted 
bv IVikhlinin et al.1 (120051) . the conversion of gas into stars 
is a crucial factor that should be taken into account. Ac- 
tually, it would be very interesting to measure whether it 
could completely explain the observed central baryon de- 
pletion on its own, or by the contrary some other physical 



mechanism (e.g. heating) must be invoked in order to ex- 
plain the observed d ensity and temperature profiles (see e.g. 
iBorgani et afll2005L for a recent discussion on this issue). 

In any case, radiative processes must obviously affect 
the physical pro perties of real cluste rs to some extent. As re- 
cently shown bv lO'Hara et al.| (I2005T) , the scaling relations of 
clusters with and without a cool core are clearly offset from 
each other, and the observed scatter can be significantly re- 
duced by introducing the peak strength (characterized by 
the central X-ray surface brightness) as an additional pa- 
rameter. 

Although such a parameter would measure a mixture 
between the intensity of cooling and the internal structure 
of the halo, a visual comparison between the simulation re- 
sults plotted in Figure|3]and the observational data shown in 
Figure|3suggests that observed systems display a somewhat 
larger scatter than our simulated clusters, which points in 
the direction that both processes may have a comparable 
contribution to the total scatter. Measurement errors, most 
notably for cluster mass estimates, also contribute to the 
scatter in the observed scaling relations, although they have 
been reported to be relatively sm all compared to the in- 
trinsic scatter JO'Hara et alj l2005). We have also neglected 
redshift evolution, which can modify the observed masses 
and luminosities by a factor H(z)/Ho, which for a ACDM 
universe amounts to about 40 per cent at z = 0.2. A rig- 
orous statistical analysis (and a larger dataset) would be 
required in order to make a quantitative assessment. 



5 CONCLUSIONS 

In this paper, the scaling relations between gas and dark 
matter mass, X-ray luminosity and emission-weighted tem- 
perature of galaxy groups and clusters have been investi- 
gated from a theoretical point of view. As a starting point, 
we have considered a relatively simple case in which the in- 
fluence of radiative processes on the observable properties of 
the ICM gas has been completely neglected. Our estimates 
of the adiabatic scaling relations h ave been comput e d from 
the polytropic models described in lAscasibar et all |2003), 
based on the results of high-resolution gasdynamical simu- 
lations. 

Our main conclusions can be summarized as follows: 

(i) Dark matter haloes are well known not to scale self- 
similarly, but according to a certain mass-concentration re- 
lation. We find that the effective polytropic index of the gas 
also varies systematically with mass, and propose the phe- 
nomenological fit 7 = 1.145 + 0.005 C200, equation (I23L to 
model the dependence of 7 on the concentration c of the 
dark matter halo. 

(ii) Given c(M) and 7(c), the whole structure of the ICM 
is fully specified by our model. It turns out that the effects of 
the varying polytropic index and concentration tend to can- 
cel out at all overdensities, yielding scaling relations that 
are well described by simple power laws whose exponents 
coincide with the self-similar prediction, and whose normal- 
izations are well fitted by a 'typical' c v i r ~ 8. 

(iii) Our model provides an excellent match to numerical 
data. The normalization of the M-T relation is significantly 
lower than previous values reported in the literature, which 
we attribute to a resolution effect. The scaling of the baryon 
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Figure 5. Observed M — Tx, F — Tx and Tx — Tx scaling relations, compared to our theoretical prediction (solid lines) and the 
one-sigma scatter expected from Ac v i r /c v j r ~ 0.3 (dotted lines). 



fraction and the L — T relation are broadly consistent with 
independent numerical work. 

(iv) Additional physics (most notably, radiative cooling 
and star formation) has an important effect on the density 
and temperature profiles of real clusters. The M — T relation 
is not severely affected, but the baryon fraction observed 
in low-mass systems is considerably below our theoretical 
prediction. This results in a lower X-ray luminosity, and it 
is ultimately responsible for the steepness of the observed 
L — T relation. On the other hand, the precise strength of 
cool cores (which cannot form in our simulations) seems to 
increase the scatter around the average scaling relations. 
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